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The experimental data from quasielastic electron scattering from are reanalyzed in terms 
of a new scaling variable suggested by the interacting relativistic Fermi gas with scalar and vector 
interactions, which is known to generate a relativistic effective mass for the interacting nucleons. 

By choosing a mean value of this relativistic effective mass = 0.8m at, we observe that most 
of the data fall inside a region around the inverse parabola-shaped universal scaling function of 
the relativistic Fermi gas. This suggests a method to select the subset of data that highlight the 
quasielastic region, about two thirds of the total 2,500 data. Regardless of the momentum and energy 
transfer, this method automatically excludes the data that are not dominated by the quasielastic 
process. The resulting band of data reflects deviations from the perfect universality, and can be 
used to characterize experimentally the quasielastic peak, despite the manifest scaling violation. 

Moreover we show that the spread of the data around the scaling function can be interpreted as 
genuine fluctuations of the effective mass M* = m*^/mN ~ 0.8±0.1. Applying the same procedure 
we transport the scaling quasielastic band into a theoretical prediction band for neutrino scattering 
cross section that is compatible with the recent measurements and slightly more accurate. 


Quasielastic electron scattering from nuclei has experi¬ 
enced a revival out of the practical need to gauge the va¬ 
lidity of the current models as apriied to neutrino scatter¬ 
ing and oscillation experiments (for recent reviews 
see refs. @-0). The conventional approach pursues a de¬ 
tailed microscopic relativistic description of the inelastic 
processes and then requires all the relevant mechanisms 
for the particular kinematics. At present there is no 
compelling model able to describe the world (e,e') ex¬ 
periments. In the case of taken as example here, 
the more than 2,500 data available spread over a huge 
(g, (jj) kinematical region, reaching well inside the rela¬ 
tivistic regime. A crucial issue is to find which electron 
data encode the maximum information to be applied to 
neutrino scattering minimizing the systematic and theo¬ 
retical uncertainties in the relevant channel (quasielastic, 
pion emission, ...). The scaling approach provides an ap¬ 
pealing and unified framework to encompass coherently 
the large diversity of data stemming from different ex¬ 
periments and kinematics. In particular the super-scaling 
approach (SuSA) has been implemented along these lines 
to predict neutrino scattering cross sections from a lon- 
gitudinal scaling function /lCV'O fitted to electron data 
[l0|. Moreover, the recent upgrade of the SuSA-v2 111 
includes nuclear effects which are theoretically-inspired 
in a particular realization of the relativistic mean field 
(RMF) theory, by an additional transverse scaling func¬ 
tion which is different from 

A peculiar feature of the RMF is that it is the only ap¬ 
proach which reproduces the experimental scaling func¬ 
tion /l(^ 0 fo'' tfi® values of q after an ad-hoc q- 
dependent shift in energy is applied [l^ . The theoret¬ 
ical origin of this phenomenological shift has not been 


well understood llj. This model incorporates a dynam¬ 
ical enhancement of lower Dirac components, which is 
transmitted to the transverse response, Rt, improving 
the agreement with experimental data. However, gauge 
invariance is violated and hence Rt still presents ambi¬ 
guities [l^. Another difficulty of the RMF and other 
finite nuclei models is that they break translational in¬ 
variance (attempts to restore it in a relativistic system 
were explored in refs. [ilEl)- 

The goal of this paper is to exploit the scaling idea 
from a novel point of view connecting the RMF with the 
universal scaling function of the relativistic Fermi gas 
(RFC) 
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Rather than constructing a yet undetermined scaling 
function we aim to propose a new scaling variable ip* 
mapping the data into a region around the above func¬ 
tion. Inspired by the fact that the mean field theory pro¬ 
vides a consistent and reasonable description of the nu¬ 
clear response in the quasielastic region (already observed 
by Rosenfelder 35 years ago [l^) for a range of kinemat¬ 
ics, we propose to start from the interacting RFG 171 
including a suitable vector and scalar potentials which 
are inferred from the data into an effective mass 
that gets reduced in the nuclear medium. The effective 
mass encodes relativistic dynamical effects relevant in 
this kinematical region, alternative to other app roaches 
like the one based on the spectral function [I1[I3. In 
fact, one of the motivations of our approach, called here 
M*-scaling (or M*S), was to provide a framework en¬ 
joying the good features of the RMF without incurring 
into the above mentioned difficulties, unvealing the m'^ 
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origin of the dynamical enhancement of both the lower 
components and the transverse response function. The 
shift of /l(V' 0 is trivially obtained as a consequence of 
the m* dependence of the quasielastic peak position. 

Of course, there have been numerous attempts to de- 

MEl, 


termine the effective mass 
details of the dynamics. Thus 


but this depends on 
by proceeding directly 


from the data we avoid specifying the mean field ex¬ 
plicitly. On the other hand a phenomenological deter¬ 
mination of m*j^ suffers from the uncertainties on the 
bulk of the data which should most signihcantly con¬ 
tribute. Therefore, we will from the beginning accept 
that this effective mass is determined up to a sensible 
uncertainty, defined precisely by a suitable selection of 
the large database, to be explained below in detail. One 
of the main advantages of this rather simple approach, 
is not only its ease of implementation, but also that we 
are free from the traditional objections regarding gauge 
invariance or PCAC violations. We expect in this way to 
account for the most relevant uncertainties regarding the 
predictive power of the model. 

We follow closely the notation introduced in Ref. (2^ . 
The quasielastic electroweak cross section is proportional 
to the hadronic tensor or response function for single¬ 
nucleon excitations transferring momentum q and energy 
w, which in the Fermi gas reads 


^ J d?p5{E'-E-u:) 


EE' 


X 2w^.n. (P'> v)0{kF - p)0{p' - kp) (2) 


where E = is the initial nucleon energy in 

the mean field. The hnal momentum of the nucleon is 
p' = p -I- q and its energy is E' = + Note 

that initial and final nucleons have the same effective 
mass m*j^. The volume V = iir^N/k^ of the system re¬ 
lated to the Fermi momentum kp and proportional to 
the number N of protons and/or neutrons participating 
in the process. Finally the electroweak interaction mech¬ 
anism is implicit in the single-nucleon tensor 

<n.(p',P) = 

ss' 


where is the electroweak current matrix element be¬ 
tween free positive energy Dirac spinors, with mass 
and normalized to uu =1. In the case of electron scat¬ 
tering we are involved with the electromagnetic current 
matrix element 


JpsiP'^P) = Us'{p') 




Q. 
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?^s(p) 

( 4 ) 

where Ff and F 2 are, respectively, the Dirac and Pauli 
electromagnetic form factors of proton or neutron. 

In the case of (e, e') the quasielastic cross section is 
written in Rosenbluth form 


da 

d^'de' 


O’Mott {vlRl + vtRt) 


( 5 ) 


where ctMott is the Mott cross section, vp = and 

vt = tan^(0/2) — Q‘^l2q‘^, with 0 the scattering angle. 
The nuclear longitudinal and transverse response func¬ 
tions are the following components of the hadronic tensor 
in a coordinate sustem with the z-axis in the q direction 
(longitudinal) 


i?L(g,w) = IF™ (6) 

RT(g,w) = IF™+IF^2 (7) 

In the RFG the nuclear response functions can be written 
in the factorized form ior K = L,T 

Rk = GKfir), (8) 

Gk = A(ZC/P + iVC/^) (9) 

Where /(t/*) is given in Eq. (1) and tp* is defined below. 
Moreover 


A = 






and the single nucleon response functions are 


Ul = — 

T 


{G*pf 


{G*pY + t{G\j) 

l+T 


Ut = 2r(Gl,f + ISAlASAi 

i + T 

where the quantity A has been introduced 
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(13) 


Dimensionless variables have been introduced measuring 
the energy and momentum in units of m'^, namely A = 
ujf2m’^, K = q/2m*pi^ t = — A^, rjp = kp/my, and 

Cf = yi + 77 I. — 1 . Note that usually 2^ these variables 
are defined with respect to the nucleon mass rriN instead 
of the m*j^. The same can be said with respect to the 
electric and magnetic form factors, that are modified in 
the medium due to the effective mass according to 



II 

(71 

1 

(14) 

rriN 


= Fi+^F2. 

(15) 


One should still stress that Fi and F 2 can depend on M* 
(^ . We stick here to the phenomenologically succsess- 
full CC2 prescription that reproduces the experimental 
superscaling function [l^. Using the CCl operator ob¬ 
tained through the Gordon reduction produces the same 
effects as in the RMF of ref. [l^ . The same modification 
of form factors in the medium was explored in ref. [ 2 ^. 
For the free form factors we use the Galster parametriza- 
tion. 
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To define the scaling variable, ip*, we first introduce the 
minimum energy allowed for a nucleon inside the nucleus 
to absorb the virtual photon (in units of m^) 


Co = Max K\l IH-A, cf — 2A 


(16) 


where ep = + r/p is the Fermi energy in units of 

The scaling variable is defined by 


r 


- r) 


(U) 


Note that ip* < 0 for A < r (the left side of the quasielas¬ 
tic peak). The meaning of ip*"^ is the following: it is the 
minimum kinetic energy of the initial nucleon divided by 
the kinetic Fermi energy. 

Starting with the experimental (e, e') cross section we 
compute the experimental scaling function /exp 


( da \ 

V dCl'de' / exp 

O’Mott {vlGl + VtGt) 


(18) 


which would correspond to the function f{ip*) in the rel¬ 
ativistic Fermi gas model. 

We summarize the results of our M*S analysis in figure 
la. We plot the experimental scaling function for the 
bulk of data 11, as a function of the scaling 
variable ip*. We take 


M* = ^ = 0.8. 
rriN 


(19) 


We see that a large fraction of the data collapse into a 
data cloud surrounding the RFG scaling function, given 
by Eq O- Other choices of m’^ are possible but the 
clustering substantially detunes from the RFG. So we in¬ 
terpret this pattern as the kinematic regions highlighting 
the effective Fermi gas behavior of data. This collapse 
of data resolves two issues simultaneously: On the the¬ 
oretical side it provides an operational definition of the 
relativistic effective mass, whereas on the experimental 
side provides an operational definition of the quasielastic 
peak behavior. 

The observed scaling is not perfect in the sense that 
the blur of data presents a finite width, but the width is 
roughly homogeneous as seen in Fig. 2. There we select 
the data that are clustered on a coarse grain scale ac¬ 
cording to a method inspired by the visual and conven¬ 
tional Gaussian low-pass filtering (Gaussian blur) (28| . 
Due to the discrete, heterogeneous and finite nature of 
the data in our case we use instead a constant weighting 
function. This function measures the density of points 
clustered above a given threshold m, inside a circle of 
radius r centered at the experimental point, plus minus 
the experimental error. In the figure we show four situ¬ 
ations corresponding to r = 0.1, and for illustration the 
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FIG. 1: Top panel (a): M* scaling analysis of the experimen¬ 
tal data of as a function of the scaling variable ip* for 
M* = m*i^/mN = 0.8 compared to the RFG parabola. Bot¬ 
tom panel (b): RFG Monte Carlo simulation of QE data with 
a Gaussian distribution of relativistic effective mass quotient 
around M* = 0.8 ± 0.1. The Fermi momentum is fixed to 
kp =225 MeV/c. 


result of applying our low pass filtering method to four 
values of TO = 20, 25,30 and 40. The parameter m mea¬ 
sures the minimum number of experimental points sur¬ 
rounding each data in the cloud. Note that we discard 
the surrounding points that don’t verify the above condi¬ 
tion. As we can see the shape defined by the data cloud, 
seen as a shaded band in the scale of the figure, presents 
a stable pattern around the relativistic Fermi gas when 
the threshold value increases, even if the number of sur¬ 
viving points decreases. This stability around the Fermi 
gas result is triggered by the chosen value of M*. Note 
that the number of data involved in these plots is around 
1,500, but the scaling violation (defined as the width of 
the shaded band) is manifest. 

This pattern in the M*S plot, that emerges as a re¬ 
alization of an universal quasielastic peak, is a global 
property of the set of data and suggests an alternative 
interpretation in terms of fluctuations of M*. One could 
propose a statistical model where each point in the cloud 
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FIG. 2: Experimental data selection in terms of the scaling 
variable tp*, obtained with different choices of the number m 
of points inside a circle with radius r = 0.1. 






samples a quasielastic event with a slightly different ef¬ 
fective mass around the mean value 0.8. This fluctua¬ 
tion does not simulate the nuclear effects beyond the im¬ 
pulse approximation (finite size effects, short-range NN- 
correlations, long-range RPA, meson-exchange currents, 
A excitation, pion emission, two-particle emission, final 
state interaction). However the fluctuations are of the 
same order of magnitude as these effects. Actually they 
are small enough to retain the points in the neighborhood 
of the quasielastic region, which could be treated pertur- 
batively in a microscopic framework beyond the RFG. 
The largest deviations of the quasielastic cloud from the 
perfect parabola occur only around its edges, where the 
Fermi gas is zero and hence the resulting signal cannot 
be accounted for by a change of 

To justify the above assumption, we carried out a cal¬ 
culation using a family of RFGs with slightly different 
m^, to generate a random point for each single exper¬ 
imental datum at the very same kinematics. Thus we 
take a random M* around the optimal mean value 0.8 in 
a Monte Carlo sampling. The results of this simulation 
are shown in Fig.l-b. To generate the pseudo data we use 
a Gaussian distribution with a width value a = 0.1, rep¬ 
resenting the fluctuation of M* = 0.8 ± 0.1, that nicely 
resembles the fluctuations seen in the cloud of the exper¬ 
imental data. This procedure automatically selects those 
pseudo data attributable to genuine quasielastic interpre¬ 
tation (based on the Fermi gas definition) and makes zero 
those kinematics that are forbidden. 


b 




FIG. 3: Total QE neutrino cross section of! per neutron 
as a function of the neutrino energy for different relativistic 
effective masses generated in a Monte Carlo simulation. The 
experimental data points are from NOMAD [J] and Mini- 
BooNE [2|. We take the axial dipole mass Ma = 1 GeV. 


Our main observation is that by choosing the optimum 
relativistic effective mass, a RFG-like scaling of the data 
can be obtained in the quasielastic region, covering more 
than 1500 data. This implies a trade off between the 
experimental uncertainty to what extent a datum is close 
to quasielastic and the importance of the physical effects 
beyond the impulse approximation that contribute to the 
quasielastic mechanism. With this procedure a way to 
estimate what information is contained in the data about 
the quasielastic peak emerges. 

From our analysis a phenomenological scaling function 
could be also obtained exactly in the same way as in the 
superscaling analysis. We have not tried to parameter¬ 
ize this function, that could be done from the data of 
fig. 2. The resulting scaling function is asymmetrical 
and very similar to the longitudinal superscaling func¬ 
tion, but with a different normalization, including the 
tail. Therefore the tail of the scaling function is a prop¬ 
erty of the quasielastic interaction. 

The information extracted here about the quasielastic 
(e, e') cross section of ^^G can be straightforwardly used 
to make readily predictions for other reactions like GG 
neutrino scattering from the same nucleus. In figure 3 
we show the calculations of the cross section as 

a function of incident neutrino energy. There we show the 
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effective RFG results with = O.Sitin- The cloud of 
points correspond to incorporating the same fluctuations 
±0.1 of M* as in the Monte Carlo simulation depicted 
in Fig. 1-b. By comparison we also show the results of 
the conventional RFG. The effective mass produces an 
enhancement of the lower Dirac components, and hence 
also of both the vector and axial transverse responses, 
and of the theoretical cross section, which thanks to the 
fluctuations becomes compatible with the data for all the 
kinematics. In our case the fluctuations of the theoretical 
band are about 10%, as naively expected from the input 
uncertainty of the effective mass. We note that the sam¬ 
pling in the lower panel of figure 3 uses a smaller binning 
of 1 MeV as opposed to the upper panel where 100 MeV 
was used instead. The clustering of these equidistant 
binnings arises naturally from the log scale. 

For high Q the vector form factors deviate from the 
conventional dipole behaviour 2^, which could affect, in 
principle, any model’s predictions. However this would 
only be appreciable in the differential cross section; the 
integrated cr, even for NOMAD kinematics, is only sen¬ 
sitive to the kinematical regions where the product of 
form factors and phase space is large. We have numer¬ 
ically checked that the contribution from above 1-2 
(GeV/c)^ is negligible, because of the rapid fall of the nu¬ 
cleon form factors. As a matter of fact one can ignore the 
electric neutron form factor completely after integration. 

Note that the set of data of unfolded energy depen¬ 
dent CCQE cross section model suffer from uncertainties 
driven by the model dependence of the neutrino energy 
reconstruction. The comparison of fig. 3 is merely in¬ 
dicative for illustration purposes of the kind of predic¬ 
tions that the present approach can provide for proper 
flux-averaged doubly differential cross sections. These 
comparison will be presented in a forthcoming publica¬ 
tion. 

At present there is no model able to reproduce the 
2,500 data points from ^^C{e, e') experiments. Due to 
the impossibility to fit the quasielastic peak or other 
regions with the experimental accuracy, in the present 
approach, we have shown that instead of making an ex¬ 
tremely detailed analysis of the particular reaction, which 
maybe well beyond the present validation possibilities, it 
is possible to isolate those data contributing to the sim¬ 
plest possible physics we are interested in, and use that 
information to make predictions with the maximum al¬ 
lowed precision, since one cannot distinguish the theoret¬ 
ical noise from the experimental signal. 
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